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I. INTRODUCTION 

During the last three decades, a large amount of work 
has been devoted to the study of nonlinear systems sub- 
ject to noise and time-dependent forces. In the course of 
these studies it has become clear that, by contrast to the 
role played by stochastic forces in linear systems, noise 
can drastically alter the response of nonlinear dynami- 
cal systems to the external forcing under some particular 
circumstances. A particularly interesting example of the 
effects of noise within the framework of signal process- 
ing by nonlinear systems is Stochastic Resonance (SR), 
i.e., the amplification of a weak input signal by the con- 
certed actions of noise and the nonlinearity of the sys- 
tem. Although discussed initially within the context of 
dynamical systems with bistable potentials the phe- 
nomenon of SR has also been found in other dynamical 
systems 0, 0, including systems of biological interest 
(see, e.g., Ref. Q and references therein). Many authors 
have studied the occurrence of SR not only in the pro- 
cessing of harmonic signals but also of aperiodic signals 
El IE SB and in the presence of both white and colored 
noise \9L lift Hl|. Furthermore, similar effects have been 
also found when a chaotic signal is used instead of noise 


The paradigmatic model in the study of SR is the over- 
damped motion of a particle in a symmetric double well 
potential driven by a periodic signal and noise. Its dy- 
namics is described by the stochastic differential equation 
(SDE) 



x(t)=-U'[x(t)]+F(t)+^(t), 



(1) 



where U'(x) is the derivative with respect to x of the 
symmetric double well potential U(x), F(t) represents 
the periodic forcing, and is a zero average Gaussian 



white noise with autocorrelation function (£(i)£(s)) = 
2DS(t— s). SR in this model can be understood in terms 
of the synchronization of two time scales: the Kramers 
time, characterizing the noise induced interwell transi- 
tions, and the time scale associated to the external driv- 
ing. Roughly speaking, when Kramers' time matches 
half the period of the external driving for a given noise 
strength D, the amplification of the weak signal F(t) is 
optimal. 

Recently, an analogous phenomenon named vibra- 
tional resonance (VR) has been shown to occur when the 
noise is replaced by a high-frequency periodic force of 
varying amplitude. Originally described by Landa and 
McClintock 0], its study has been also addressed by 
other authors and, from different points of view, in ex- 
citable 0] , spatially extended , and bistable systems 
[wL [ItJ ■ In Ref. [13 , a brief numerical study of the ef- 
fects of additive noise on VR has been presented, but a 
more detailed investigation of this topic is still lacking. In 
this paper, we undertake an analytical study of the prob- 
lem that will provide explicit expressions for the relevant 
observables. Our analytical results are compared with 
numerical simulations performed on the model described 
in Ref. 0. 



II. DESCRIPTION OF THE MODEL AND 
CHARACTERIZATION OF VR 

We consider a system described by the SDE (|TJ, with 
U(x) being the symmetric quartic potential in dimen- 
sionless form 

4 2 

U(x) = ^-^, (2) 

and F(t) the bichromatic force 

F(t) = Acos(fit) + Nnrcos(Nnt + ip). (3) 
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Here, for reasons that will become clear in the next sec- 
tion, we have introduced in the second monochromatic 
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force an arbitrary initial dephasing, (p. The parameter N 
is chosen to be a positive integer, so that F(t) is periodic 
with the same period T = 2-k/VL as the first monochro- 
matic force. By r we have denoted the ratio of the ampli- 
tude of the second force to its frequency. This ratio r is 
assumed to be of the same order as the parameters char- 
acterizing the potential U(x). We are interested in situa- 
tions in which the parameters N fir and N appearing 
in the second monochromatic force are much larger than 
the rest of the parameters in the problem. In this sense, 
we will say that N£lr cos(N fit + <p) is a strong, high- 
frequency monochromatic force. This situation can be 
formally achieved by taking the limit N — > oo, with the 
ratio r kept fixed. 

The corresponding Fokker-Planck equation (FPE) for 
the probability density P(x,t) reads 



-P(x t) = — 

at 1 ' ' dx 



D^- + U'{x)-F{t) 



P(x,t). (4) 



The analysis of this equation is simplified by making use 
of two important theorems: the H-theorem, which en- 
sures the existence of a unique long time distribution 
function P oc (x,t) 0, 0, and the Floquet theorem, 
which guarantees that Poc(x,t) is periodic in time with 
the same period, T, as the external force [2(j. Hence- 
forth, we will assume that the relaxation transient stage 
has ended and, consequently, the long time distribution, 
P oc (x,t), has been reached. 

The first moment of the probability distribution can be 
used to characterize the system response to the external 
driving. As a consequence of the above mentioned theo- 
rems, its long time limit, (^(t))^, is a periodic function 
of time with period T and, therefore, it can be expanded 
in Fourier series as follows 



This parameter represents the ratio of the amplitude as- 
sociated to the first harmonic to A, and it is directly re- 
lated to the spectral amplification r\ through the expres- 
sion Q = ^ff\ (see Ref. |2|). Notice that the above defi- 
nition of the parameter Q differs by a factor 2 from the 
definition used in the theoretical approach in Ref. |l7j . 
and by a factor A~ 1 from the one used in the numerical 
results appearing in the same reference. In the determin- 
istic case, D = 0, it has been shown that Q goes through 
a maximum as the intensity of the high-frequency force 
is increased. This maximum is approximately localized 
around the value r = •v/2/3 of the ratio of the ampli- 
tude of the second force to its frequency 0, 0] . This 
behavior resembles that of SR, with the noise strength 
D playing the role of the intensity of the high-frequency 
force. 

When white noise is added, the situation is different. 
More precisely, the two main effects of increasing the 
value of the noise strength D are that the value at the 
maximum of the resonance curve decreases and, at the 
same time, its location is shifted towards lower values 
of the high-frequency amplitude. Even more, for large 
enough values of D, the VR phenomenon completely dis- 
appears 0. These effects can be qualitatively under- 
stood as a consequence of the fact that the white noise 
provides an input to the system with contributions to 
all the frequencies. On the one hand, the fraction corre- 
sponding to the high-frequency region advances the ap- 
pearance of the maximum. On the other hand, the re- 
maining contribution masks the high-frequency force, de- 
creasing its relative importance . To shed some more 
light on this question, an extension of the theoretical ap- 
proach put forward in Ref. |17| to the case in which white 
noise is present would be desirable. This is the main aim 
of the next section. 
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f(t)>oo = X " 
n— — oo 

oo 

= X + 2j2\Xn\cos{nflt -<£„), (5) 



where 



1 f T 

Xn = T J dt ( X ^< 



and the phases 4> n have been chosen so that 
cos (f>n — Re I 
sin <b„ = — Im 



X n \ 
X 



(6) 

(7) 
(8) 



The observable that has been usually considered to 
quantify the VR phenomenon is the parameter 



Q = 2 



A 



(9) 



III. THEORETICAL APPROACH 

In this section, we will obtain analytical expressions for 
the parameter Q based on three approximations of differ- 
ent nature. The first two ones are simply generalizations 
to the noisy problem of those carried out in the theoreti- 
cal approach developed in the absence of noise • The 
third approximation is specific to the noisy problem. For 
clarity in the exposition, as well as in the discussion of 
their validity conditions, we will present them separately. 



A. Derivation of the effective dynamics 

As a consequence of Eqs. © and 10} , the time deriva- 
tive of the probability density P(x, t) diverges as N in the 
limit N — > oo. To avoid this divergence, it is convenient 
to extract the fast dependence from x(t), and define the 
new stochastic process 



y(t) = x(t) -rsin(Nnt + ip). 



(10) 
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The Fourier components of (x(f))oo and (y(f))oo are re- 
lated according to 



X n + '-k (e 



-N 



(11) 



where Y n is defined from Eq. by replacing (a;(f))oo 
by (y(t))oo- Therefore, as iV 7^ ±1, we can substitute the 
coefficient X\ by Y\ in the definition of the parameter Q 
[see Eq. ©]. 

The advantage of using the process y(t) instead of x(t) 
becomes clear after writing the FPE for its probability 
density, which reads 



d_ 

at 



d_ 

dy 



D—+U'(y,t,<p) 
dy 



- Acos(ttt) 



r(y,t^). 



(12) 



Here, we have introduced the new potential U(y, f, tp) := 
U[y + r sin (TV Slt + ip)], and the dependence of the func- 
tions on the phase tp has been written explicitely. From 
Eq. I|12|) and the definition of U (y, t, tp) it follows that 
the time derivative of V(y,t,tp) is at most of order 1 
as N — > 00. Therefore, a large number of oscillations 
of the function r sin(7V Q t + ip) appearing in the defini- 
tion of lA{y 1 f, ip) take place before a significant change in 
P(y,t,ip) occurs. As a consequence, it is to be expected 
that, for N ^> 1, V(y,t,tp) is almost independent of the 
phase tp. If we define the phase average of an arbitrary 
function f(y,t,tp) as 
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f{y,t) = f{y,t 1 ip):=— dipf(y,t,tp), (13) 
27T Jo 

then, for N ^> 1 and any value of ip G [0, 2tt], one can ap- 
proximate V(y,t,tp) ps V(y,t). Furthermore, the decou- 
pling approximation f(y,t,tp)V(y,t,tp) w f(y,t) V(y,t) 
also holds. Carrying out the phase average in Eq. H12|) 
and using the decoupling approximation, one obtains 



d_ 

m 



V{y,t) 



d_ 

dy 



D 



d_ 

dy 



Kf f (y) 



Acos(Clt) 



V(y,t), 

where we have in troduced the effective potential 
U e ff(y) '.— U(y,t, tp), which is clearly time independent. 
This averaging procedure has been previously used in the 
study of thermal activation in bistable systems under the 
influence of a periodic force with moderate to large fre- 
quency |2lj |. It has also been used in the explanation 
of the enhancement observed in the tunnel splitting of 
a quantum bistable system when a very high frequency 
driving is applied 22] . 

It also follows from the above considerations that, for 
N Z>1, (v(t))oo « (y(t))W f) := f^dyyV^t), 
where Voc(y, t) is the long time distribution correspond- 



ing to the FPE i|14|) . In particular, within this approxi- 
mation, the parameter Q is given by 



Q 



(e//)| 



.4 



(15) 



where is obtained from Eq. JBJ) with n = 1, by 

replacing (x(t))oo by (y(*))£> ■ 

The explicit calculation of U e f / (y) for the potential in 
Eq. J3J) leads to 



u eff(y) = \~ a ( r )y' 



(16) 



where we have dropped an irrelevant constant and have 
introduced the quantity air) = 1 — 3r 2 /2. From this re- 
sult, it is clear that the stability of this effective p oten tial 
depends on the ratio r. More precisely, if r < y/2/3 the 
potential is bistable, whereas if r > a/2/3 it is monos- 
table. Thus, an increase in r leads to a decrease in the 
effective barrier height and, eventually, to its disappear- 
ance. 

In summary, we have shown that, for N ^> 1, the 
time evolution of the original stochastic process x(t) can 
be approximately described in terms of the dynamics 
of a Brownian particle moving in the effective potential 
U e ff(y) and in the presence of noise and the monochro- 
matic force Acos(Qt). 



B. Linear Response Theory 

Henceforth, following the approach in Ref. 0]> we will 
assume that the amplitude A is small enough, so that 
Linear Response Theory (LRT) provides a good descrip- 
tion of the effective dynamics obtained in the previous 
subsection. An extensive study of the validity conditions 
of LRT can be found in Refs. |H Within LRT 
the parameter Q in Eq. I|15(l is given by [2j 



where 



Xeffity 



\xef f m 



-iSlt 



(17) 



(18) 



is the value at to = Q of the Fourier transform of the 
response function, x e ff(t), corresponding to an over- 
damped Brownian particle moving in the effective po- 
tential U e ff(y). This response function obeys the well- 
known Fluctuation-Dissipation Theorem (FDT) [26l I27I 

M 



Xeff(t) = -^±K eff (t), 



(19) 



where 77(f) is the Heaviside step function, and K e ff(t) = 

(y{t)y(Q))e(/^ is the equilibrium autocorrelation func- 
tion of the effective system described by the FPE 114|) 
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in the absence of external driving. Thus, the stochastic 
process y(t) appearing in the definition of K e tf{t) is a 
solution of the SDE 



m = -u' eff [ y (t)}+m- 

Use of Eq. flU in Eq. (JTSJ yields 



(20) 



dtK eff (t) sin (nt) 



KeffjO) 



(21) 



lm[xeff(m = --J dt K eff (t) cos (nt). (22) 

These two equations, together with Eq. I|17|) . allow us to 
express Q( LRT ^ in terms of the equilibrium autocorrela- 
tion function K e ff(t). 

For the nonlinear potential in Eq. 1)16(1 , explicit expres- 
sions for K e ff(t) are unknown, so that this correlation 
function must be evaluated either numerically or by re- 
sorting to suitable approximations. Before applying to 
our problem approximate techniques discussed in the lit- 
erature, it is convenient to reduce the SDE (|20|l to a more 
standard form. In order to do so, let us rescale the co- 
ordinate and time, and define the new stochastic process 
y(t) = KOI-^da^j-if) (for r ^ y^2/3). Then, 
from Eq. 1121 III , it is easy to verify that y (i) fulfills the 
SDE 



where we have introduced the rescaled potential 



u(v) 



fj 4 



3gn[o(r)] 



y 2 



2 ' 



(23) 



(24) 



and the rescaled Gaussian white noise with zero 
average and autocorrelation function (£m£(t')) = 
2 |a(r)|~ 2 Dd(t — t'). From these considerations, it is 
straightforward to prove that the equilibrium autocorre- 
lation function of the original process, K e ff(t,D), and 

that of the rescaled one, K (t, D), are related by 



K eff (t,D) = \a(r)\K (\a(r)\t, 



D 



\a(r)\< 



(25) 



Analogously, Eqs. 0, HU and lead to the fol- 
lowing relation between the Fourier transform of the re- 
sponse function of the original process, | Xef f (^ 5 D) | , and 
that of the rescaled one, \x(Q, D) | , 



\x eff (n,D)\ = 



\a(r)\ 



D 



\a(r)\> \a(r)\* 



(26) 



C. Weak Noise Approximation 

The expressions (|25|l and l|26|) allow us to evaluate 
K e ff(t,D) and |x e //(S7, D)| from the rescaled functions 

K(t, D) and |x(0,D)|. The asymptotic behaviors of 
these rescaled functions for small values of the noise 
strength D have been widely studied in the literature. 
To summarize the results, we will consider separately two 
cases. 

Firstly, if r < y^2/3, then U{y) = j/ 4 /4 - y 2 /2 is the 
archetypal symmetric quartic double-well potential ex- 
pressed in dimensionless form. For this bistable poten- 
tial, we can use the two- mode approximation |29j . This 
approximation is based on the existence of a clear-cut 
separation between the time scales associated to inter- 
well and intra-well motions, and it is expected to be valid 
for small values of the noise strength D. Within this ap- 
proximation, K(t,D) and |x(f2, .D)| are given by 

K(t, D) = (ft (D) exp [ - Ax (D) t\ + g 2 (D) exp [ - at] , 

(27) 

and 



\x(&,D)\ = 



1 f 2ag 1 (D)g 2 (D)\ 1 (D)[a\ 1 (D) + Q 2 ] 



D 



[x 2 (D) + n 2 ] (a 2 + n 2 ) 



+ 



9l{D)\\(D) , g 2 (D)a 2 



1/2 



X 2 (D)+n 2 a 2 +0 2 



(28) 



In the above expressions, the parameter a — U"(±l) = 2 
is the curvature at the minima of the potential U(y), and 
Ai(D) is the the smallest non- vanishing eigenvalue of the 
Fokker-Planck operator corresponding to the SDE J2HJ- 
In the steepest-descent approximation, this eigenvalue is 
given to leading order in D by 



Ai D) 



V2 



exp 



1 

41) 



The weights g\ (D) and g 2 (D) read 



92(D) 



[Xi(D) +l]m 2 (D) -m 4 (£>) 



a 



gi (D) = m 2 (D) - g 2 (D) 



(29) 



(30) 



(31) 



where 
m 2n (D) 



-u(y)/D 



f^dye- 11 ^)/ 

/2 r(n + l/2)D_ n _ 1/2 [- [2D) 1/2 ] 



(2D) 



T(l/2)D_ 1/2 [-(2D) 1/2 ] 



(32) 
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are the nonvanishing moments of the equilibrium distri- 
bution. Using the asymptotic expansion of the parabolic 
cylinder functions for small values of D, these weights 
can be approximated by gi(D) ~ 1 — (a + l)D/a and 
92(D) w D/a to first order in D. 

Secondly, if r > y^2/3, then U(y) = y 4 /4 + y 2 /2 is a 
monostable potential. In this case, for small values of D, 
the autocorrelation function K(t, D) is given by j^lj 



K(t,D) = D cxp [- |Ai(Z3)|t], 



with 



Ai(5) 



1 - 6D 



1/2 



(33) 



(34) 



Consequently, the Fourier transform of the response func- 
tion reads 



\x(^D)\ = 



[Af (D) + fi 2 



1/2- 



(35) 



To evaluate K e f f(t, D) and |x e j/(f2, D)\ one simply in- 
troduces either the expressions (|27|) and l|28[l or (|33() and 
(|33|l into Eqs. (|2"5|l and lj2U|) . Finally, within this approx- 
imation, the parameter Q will be evaluated by replac- 
ing in Eq. I|17|l the result obtained for \xeff(^, D) | . It 
is important to emphasize that the approximate expres- 
sions lj2*7|) . l|33|> and have been obtained under 
the assumption that D is a small parameter. Therefore, 
the resulting expressions for K e ff(t,D), |x e //(f2, D)| 
and Q are expected to be valid only for small values 
of D\a(r)\~ 2 . Taking into account that, for D ^ 0, 
the parameter D\a(r)\~ 2 diverges at r = ^/2/3, then 
it is clear tha t, ev en for D <C 1, there exists a region 
around r = in which the resulting expressions for 
K e ff(t,D), \x e ff(il,D)\ and Q are not applicable. In 
particular, the expression for Q is expected to provide a 
successful explanation of the VR phenomenon whenever 
its maximum takes place at a value of r outside the above 
mentioned region. 

From the above considerations, it is also possible to ob- 
tain an expression for the location of the maximum char- 
acterizing the VR phenomenon in the presence of white 
noise. According to the numerical results put forward in 
Ref. |17| . this maximum is located at a value ru such 
that Tm < a/2/3. Then, to evaluate the parameter Q in 
the neighborhood of tm , one has to consider Eq. Q28JI and 
follow the procedure described in the above paragraph. 
As the expression thus obtained is rather cumbersome, it 
is convenient to introduce a further approximation. More 
precisely, we will assume that, around tm, the contribu- 
tion due to the intrawell motion is negligible, so that just 
the second term in Eq. (|28|l with gi(D) = 1 is present. 
Therefore, the parameter Q as a function of r obtained 
from Eqs. (|17|l and Ij26(l is given by 



Q(r) = 



a(r) 



Ai [D/a 2 (r)} 



D 



{\ 2 [D/a 2 (r)} + [n/a(r)} 2 } 



1/2 ' 



(36) 



for r < ^2/3. The location of the maximum, cm, is 
readily determined from the equation Q'{tm) = 0, where 
Q'{r) represents the derivative of Q(r) with respect to r. 
Consequently, from Eq. I|36|l it follows that tm obeys the 
transcendental equation 



A? 



D 



_a 2 {r M ) 



= 2 



a(r M ) 



a 2 {rM) 
AD 



1 



(37) 



After making the replacements Z?/a 2 (rjv/) — > D and 
Q/a(rM) — > H in the above equation, the resulting ex- 
pression resembles, except for a factor of 2, the transcen- 
dental equation obtained within the two-state model to 
determine the noise intensity that maximizes the spectral 
amplification (see, e.g., Ref. and references therein). 
The difference in the factor of 2 is due to the fact that, 
in that case, only the noise intensity is varied, whereas in 
our case, by varying r, we are simultaneously modifying 
the frequency and the noise intensity. 

The transcendental equation (|37|l provides an ex- 
planation for the numerically observed behavior of 
the location of the maximum as a function of D 
(see Sec. |HJ|. For a fixed value of D, the func- 
tion \\ [D/a 2 {r)\ is a monotonously increasing func- 
tion of r in the interval [0, -y/2/3). Also, for 
fixed values of D and Q, 2 [n/a{r)f [a 2 (r)/(4L>) - l] 
is a monotonously decreasing function of r in the 
same interval. Thus, the existence of a solution of 
Eq. J39 in the interval (0, y/2~j?,) requires the inequality 
\\ [D/a 2 (0)] = \\ (D) < 2 [fi/a(0)] 2 ["W^) - l] = 
2 S7 2 [1/(41?) — 1]. Therefore, the same inequality is re- 
quired for the existence of VR. If we now keep fixed r 
and f2 and look at both sides of Eq. (|3Tjl as a func- 
tion of D, we notice that 2 p,/a{r)] 2 [a 2 (r)/(4£>) - l] 
is a monotonously decreasing function of D, whereas 
A 2 [-D/a 2 (r)] is a monotonously increasing function of 
D. We then conclude that the abscissa at the intersec- 
tion point between these two functions, tm, is shifted 
towards the left as D increases. This is exactly the be- 
havior observed numerically in Ref. 0. Obviously, VR 
disappears for values of D greater than the critical value 
D c for which tm = 0- This critical value fulfills the tran- 
scendental equation 



\\ [D c ) = 2fi 2 



1 

4D~c 



1 



(38) 



IV. 



COMPARISON WITH THE NUMERICAL 
SOLUTION 



In this section, we compare the results obtained for 
the parameter Q by means of the numerical solution of 
the SDE 0J, with the analytical expressions obtained in 
the previous section. In order to evaluate numerically 
the parameter Q, a large enough number of stochastic 
trajectories, x^(t), have been generated by integrating 
the SDE for every realization of the white noise £(£), 



6 



starting from a given initial condition xq . The numerical 
solution of the SDE (JIJ has been carried out by using the 
algorithm put forward in Ref. [^Il (see also the Appendix 
in Ref. |25|). After allowing for a relaxation transient 
stage, we evaluate the long time average by 



(x(t)) c 



1 M 

-Y 



(*), 



(39) 



k=l 



where M is the number of stochastic trajectories consid- 
ered. Then, the parameter Q is obtained from Eqs. © 
and JfJJl, with n = 1, by numerical quadrature. 

In Fig. ^ we depict the dependence of Q vs. r. The 
following set of parameter values has been considered: 
N = 50, ip = 0, n = 3.4 x 1CT 2 , A = 1.1 x 1(T 2 , 
D = 1.34 x 1(T 2 in panel (a), D = 6.72 x 1CT 2 in panel 
(b), D = 1.34 x KT 1 in panel (c), and D = 4.03 x 10" 1 in 
panel (d). They correspond to four representative cases 
appearing in Fig. 6 in Ref. ^tJ- Notice that the param- 
eter values appearing in that reference must be rescaled 
in order to reduce the bistable potential used there to 
the standard dimesionless form in Eq. (J2J). With circles 
we have plotted the results obtained from the numerical 
solution of the SDE JIJ together with Eqs. © and 10. 
With solid lines we have depicted the analytical values of 
Q provided by Eqs. J2| and either Eqs. 

for r < y/2/3, OT E q s - <E3> and (Q), for r > y/2/3. The 
location of the critical value r — y/2/3 is indicated by 
the vertical dotted lines. In addition, in panels (a) and 
(b) , we have depicted with vertical dashed lines the loca- 
tions of the maxima predicted by Eq. I|37l) . The critical 
value of the noise strength at which the VR phenomenon 
disappears, obtained from Eq. I|38[) . is D c = 1.154 x 10 _1 . 
The noise strength values in panels (c) and (d) are greater 
than D c , so it is to be expected that the VR phenomenon 
is not present in these cases, as it is confirmed numeri- 
cally. 

A glance to Fig. 1 reveals that the analytical expres- 
sions provide a good description of the main features of 
the VR phenomenon in the presence of white noise, at 
least if the ratio r is far enough away from the criti- 
cal value r — y/2/3. In particular, the analytical re- 
sults describe correctly the shift of the location of the 
maximum towards lower values of r, as well as the de- 
crease of the value of Q at the maximum, as the noise 
strength increases. It is important to emphasize that, 
as we mentioned in Sec. IIII CI the analytical results are 
not applicable in a neighborhood of the critical value 
r = 1/2/3, as a consequence of the divergence of the 
effective noise strength D/[a(r)] 2 . The appearance of 
the vertical asymptote at this critical value makes this 
failure evident. Notice that, in panel (a), the agreement 
between the numerical and analytical results around the 
maximum gets worse. This is due to the proximity of the 
location of this maximum to the critical value r — 
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FIG. 1: Dependence of the parameter Q vs. r corresponding 
to the noise strength values D = 1.34 x 10 -2 [panel (a)], 
D = 6.72 x 10" 2 [panel (b)], D = 1.34 x KT 1 [panel (c)], 
and D = 4.03 x 10 _1 [panel (d)]. The rest of the parameter 
values are: V = 50, ip = 0, fi = 3.4 x 10~ 2 , A = 1.1 x 10 -2 . 
With circles we have plotted the results obtained from the 
numerical solution of the SDE Q together with Eqs. (5J and 
© . With solid lines we have depicted the analytical values of 
Q provided by Eqs. JT7J, © and either Eqs. for 
r < y/2/3, or Eqs. ^3 and JSU, for r > V / 2A The dotted 
lines indicate the location of the critical value r — ^/2/3, and 
the dashed lines the location of the maximum predicted by 
Eq. (EZJ. 



V. CONCLUSIONS 

We have studied the motion of a particle in a bistable 
potential in the presence of white noise and an exter- 
nal bichromatic force, within the context of the VR phe- 
nomenon. Analytical expressions for the parameter Q 
which quantifies this phenomenon have been obtained 
based on three simplifying approximations of different 
nature. Firstly, the exact dynamics of the system has 
been reduced to that of a Brownian particle moving in 
an effective potential and under the influence of white 
noise and one of the monochromatic forces. In order 
to make this approximation, we have assumed that the 
other monochromatic force is a strong, high-frequency 
field. Secondly, we have applied LRT to express the pa- 
rameter Q in terms of the equilibrium autocorrelation 
function associated to the effective dynamics in the ab- 
sence of the external driving. Finally, we have obtained 
analytical expressions for Q within the weak noise ap- 
proximation. To do so, it has been necessary to consider 
two disjoint regions of values of the ratio r separated 
by the critical value r = y^2/3. We have also obtained 
a transcendental equation for the location of the maxi- 
mum of Q(r) which characterizes the VR phenomenon, as 



7 



well as for the condition under which this phenomenon 
dissapears. Comparison with the numerical solution of 
the original SDE shows that these analytical expressions 
provide a good quantitative description of the VR phe- 
nomenon, at least if it takes place outside a critical region 
around r = -y/2/3. With this work, we have attempted to 
provide a theoretical complement to the numerical study 
of the influence of additive noise on VR carried out in 
Ref. [13 • 
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